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Abstract 

We present an analysis of the global neutrino oscillation data in terms of four- 
neutrino mass schemes. We find that the strong preference of oscillations into active 
neutrinos implied by solar+KamLAND as well as atmospheric neutrino data allows 
to rule out (2+2) mass schemes, whereas (3+1) schemes are strongly disfavoured 
by short-baseline experiments. In addition, we perform an analysis using recent 
data from cosmology, including CMB data from WMAP and data from 2dFGRS 
large scale structure surveys. These data lead to further restrictions of the allowed 
regions for the (3+1) mass scheme. 



1 Introduction 

The neutrino oscillation interpretations of the solar 0, 0] and KamLAND_P] neutrino ex- 
periments, atmospheric 0, 0| neutrino data, and the LSND experiment p] require three 
neutrino mass-squared differences of different orders of magnitude pj. Since it is not pos- 
sible to obtain this within the Standard Model framework of three active neutrinos it has 
been proposed to introduce a light sterile neutrino to reconcile all the experimental 
hints for neutrino oscillations. Here we present an analysis of the global neutrino oscilla- 
tion data in terms of four-neutrino mass schemes, including data from solar, KamLAND 
and atmospheric neutrino experiments, the LSND experiment, as well as data from short- 
baseline (SBL) experiments 0, HI U| and reactor experiments [f3] reporting no evidence 



for oscillations. We find that for all possible types of four-neutrino schemes different 
sub-sets of the data are in serious disagreement and hence, four-neutrino oscillations do 
not provide a satisfactory description of the global oscillation data including LSND. The 
details of our calculations can be found in Refs. 0, 14 , ^ . 
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Figure 1: The six four-neutrino mass spectra, divided into the classes (3+1) and (2+2). 

2 Notations and approximations 

Four-neutrino mass schemes are usually divided into the two classes (3+1) and (2+2), as 
illustrated in Fig. [T] We note that (3+1) mass spectra include the three-active neutrino 
scenario as limiting case. In this case solar and atmospheric neutrino oscillations are ex- 
plained by active neutrino oscillations, with mass-squared differences Amg OL and Am| TM , 
and the fourth neutrino state gets completely decoupled. We will refer to such limiting 
scenario as (3+0). In contrast, the (2+2) spectrum is intrinsically different, as there 
must be a significant contribution of the sterile neutrino either in solar or in atmospheric 
neutrino oscillations or in both. 

Neglecting CP violation, in general neutrino oscillations in four-neutrino schemes are 
described by 9 parameters: 3 mass-squared differences and 6 mixing angles in the unitary 
lepton mixing matrix. Here we use a parameterisation introduced in Ref. [l3||, which 
is based on physically relevant quantities: the 6 parameters Am^, 6 S0L , Am^ TM , # ATM , 
^ m LSND' ^lsnd are similar to the two- neutrino mass-squared differences and mixing angles 
and are directly related to the oscillations in solar, atmospheric and the LSND experi- 
ments. For the remaining 3 parameters we use rj s , r\ e and d^. Here, r) s (rj e ) is the fraction of 
v s iye) participating in solar oscillations, and (1 — d^) is the fraction of participating in 
oscillations with Am| TM (for exact definitions see Ref. 13])- For the analysis we adopt the 
following approximations: (1) We make use of the hierarchy Am^ <C AmJ TM <C AmJ SND . 
This means that for each data set we consider only one mass-squared difference, the other 
two are set either to zero or to infinity. (2) In the analyses of solar and atmospheric data 
(but not for SBL data) we set rj e = 1, which is justified because of strong constraints from 
reactor experiments^, 12l |. 

Due to these approximations the parameter structure of the four-neutrino analysis 
gets rather simple. The parameter dependence of the four data sets solar, atmospheric, 
LSND and NEV is illustrated in Fig. |21 The NEV data set contains the experiments 
KARMENlj, CDHSQ, Bugey[II| and CHOOZ/Palo Verde 0, reporting no evidence 
for oscillations. We see that only rj s links solar and atmospheric data and d M links atmo- 
spheric and NEV data. LSND and NEV data are coupled by AmJ SND and 6 LSND . 
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3 (2+2): ruled out by solar and atmospheric data 

The strong preference of oscillations into active neutrinos in solar and atmospheric os- 
cillations leads to a direct conflict in (2+2) oscillation schemes. We will now show that 
thanks to the latest solar neutrino data (especially from SNO0]) in combination with the 
KamLAND experiment [3] , and the improved SK statistic on atmospheric neutrinosj^] the 
tension in the data has become so strong that (2+2) oscillation schemes are essentially 
ruled out. 1 

The thin lines in the left panel of Fig. El show the A\ 2 of latest solar data as a function 
of rj s , the parameter describing the fraction of the sterile neutrino participating in solar 
neutrino oscillations. The 99% CL bounds r/ s < 0.44, when the SSM constraint on the 
8 B-flux is included, and r) s < 0.61 for free 8 B-flux reflect the strong preference for active 
neutrino oscillations of solar data. Recently, the outstanding results of the KamLAND 
reactor experiment [3] confirmed the LMA solution of the solar neutrino problem (lR Eil] . 
Apart from this very important result the sensitivity of KamLAND to an admixture of a 
sterile neutrino is rather limited (see think lines in Fig. |HJ). The combined analysis leads 
to the 99% CL bound i] s < 0.5 for 8 B free, whereas in the SSM constrained case the bound 
is unchanged. 

In contrast, in (2+2) schemes atmospheric data imply r} 3 > 0.65 at 99% CL, in clear 
disagreement with the bound from solar data. In the right panel of Fig. |3] we show 
the x 2 f° r solar data and for atmospheric combined with SBL data as a function of rj s . 
Furthermore, we show the x 2 °f the global data defined by 

X 2 (v s ) = AxloM + Ax 2 ATM+SBL (Vs) ■ (1) 

From the figure we find that only if we take both data sets at the 99.95% CL a value 
of rj s exists, which is contained in the allowed regions of both sets. This follows from 
the % 2 -value Xpc = 12.2 shown in the figure. In Refs. 
statistical method to evaluate the disagreement of different data sets in global analyses. 
The parameter goodness of fit (PG) makes use of the x 2 defined in Eq. (jlj). This criterion 
evaluates the GOF of the combination of data sets, without being diluted by a larg e 
number of data points, as it happens for the usual GOF criterion (for details see Ref. jl9j|). 
We find Xpg = Xmm = 23.5, leading to the marginal PG of 1.3 x 10~ 6 . We conclude that 

1 Details of our analyses of the solar, KamLAND and atmospheric neutrino data can be found in 
Refs. moi. For an earlier four-neutrino analysis of solar and atmospheric data see Ref. . 



14, [ltl we have proposed a 
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(2+2) oscillation schemes are highly disfavoured by the disagreement between the latest 
solar and atmospheric neutrino data. This is a very robust result, independent of whether 
LSND is confirmed or disproved. 2 



4 (3+1): strongly disfavoured by SBL data 



It is known for a long time[2l|, |22j that (3+1) mass schemes are disfavoured by the 
comparison of SBL disappearance data0, 11] with the LSND result. In Ref. |23| we have 
calculated an upper bound on the LSND oscillation amplitude sin 2 26* LSND resulting from 
NEV and atmospheric neutrino data. From Fig. 0] we see that this bound is incompatible 
with the signal observed in LSND at the 95% CL. Only marginal overlap regions exist 
between the bound and global LSND data if both are taken at 99% CL. An analysis in 
terms of the parameter goodness of fit shows that for most values of Am 2 SND NEV and 
atmospheric data are compatible with LSND only at more than 3a, with one exception 
around Am 2 SND ~ 6 eV 2 , where the PG reaches 1%. These results show that (3+1) 
schemes are strongly disfavoured by SBL disappearance data. 



5 Comparing (3+1), (2+2) and (3+0) hypotheses 

With the methods developed in Ref. we are able to perform a global fit to the oscilla- 
tion data in the four-neutrino framework. This approach allows to statistically compare 
the different hypotheses. Let us first evaluate the GOF of (3+1) and (2+2) spectra with 
the help of the PG method described in Ref. [lij . We divide the global oscillation data 

2 Sub- leading effects beyond the approximations adopted here should not affect this result significantly. 
Allowing for additional parameters to vary might change the ratio of some observables |2d) . however, we 
expect that the absolute number of events relevant for the fit will not change substantially. 
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Figure 4: Upper bound on sin 2 2# LSND from NEV and atmospheric neutrino data in (3+1) 
schemes |23| compared to the allowed region from global LSND data0 and decay-at-rest 
(DAR) LSND data0. 





SOL ATM LSND NEV 


Xpg 


PG 


(3+1) 
(2+2) 


0.0 0.4 7.2 7.0 
14.8 6.7 2.2 9.7 


14.6 
32.4 


5.6 x 10~ 3 
1.6 x 10~ 6 



Table 1: Parameter GOF and the contributions of different data sets to Xpg m (3+1) and 
(2+2) neutrino mass schemes. 



in the 4 data sets SOL, ATM, LSND and NEV. Following Ref. (lj we consider 
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where A Xx = Xx ~ (Xx)mm {X = SOL, ATM, NEV, LSND). In Tab. □ we show the 
contributions of the 4 data sets to Xpg = Xmin f° r (3+1) and (2+2) oscillation schemes. 
As expected we observe that in (3+1) schemes the main contribution comes from SBL data 
due to the tension between LSND and NEV data in these schemes. For (2+2) oscillation 
schemes a large part of Xpg comes from solar and atmospheric data, however, also SBL 
data contributes significantly. This comes mainly from the tension between LSND and 
KARMEN |24l] . which does not depend on the mass scheme and, hence, also contributes 
in the case of (2+2). Therefore, the values of Xpg * n Tab. Q for (2+2) schemes are higher 
than the one given in Sec. 01 where the tension in SBL data is not included. 

The parameter goodness of fit is now obtained by evaluating Xpg f° r 4 DOFflij. 
This number of degrees of freedom corresponds to the 4 parameters r? s , d^, 9 LSND , Am 2 SND 
describing the coupling of the different data sets (see Eq. (j2J) and Fig. |2J). The best GOF 
is obtained in the (3+1) case. However, even in this best case the PG is only 0.56%. The 
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PG of 1.6 x 10~ 6 for (2+2) schemes shows that these mass schemes are essentially ruled 
out by the disagreement between the individual data sets. 

Although we have seen that none of the four-neutrino mass schemes can provide a 
reasonable good fit to the global oscillation data including LSND, it might be interesting 
to consider the relative status of the three hypotheses (3+1), (2+2) and the three-active 
neutrino scenario (3+0). This can be done by comparing the x 2 values of the best fit 
point (which is in the (3+1) scheme) with the one corresponding to (2+2) and (3+0). 
First we observe that (2+2) schemes are strongly disfavoured with respect to (3+1) with 
a Ax 2 = 17.8, implying a CL of 99.87% for 4 DOF. Further we find that (2+2) is only 
slightly better than (3+0), which is disfavoured with a A% 2 = 20.0 with respect to (3+1). 



6 Adding information from cosmology 

Recently a precision measurements of the cosmic microwave background (CMB) radiation 
has been published by the WMAP collaboration |25j. Under certain assumptions about 
the cosmological model and in combination with other cosmological data a stringent upper 
bound on the sum of the neutrino masses ^Zm v < 0.7 eV is obtained 25l |26j. Since this 



bound is of the same order as ,/AmL it will lead to further restrictions in four- neutrino 



LSND 

mass schemes |2 71] . To combine this information with data from oscillation experiments we 
use the results of an analysis performed by S. Hannestad[28|. There, WMAP|25T] results 
have been combined with other CMB data) 29], data from the large scale structure su rvey 
2dFGRS[30j, the HST Hubble key project [3jJ, and Type la supernovae observations [321] . 
It has been found that a non-zero m v can be compensated partially by an increase in the 
effective number of neutrino species. Therefore, the original bound obtained in Ref . 25|| , 
which was calculated assuming three neutrinos, cannot be applied in a straightforward 
way. In the case of four neutrinos a somewhat weaker bound of J2 m u < 1-38 eV (95% 
CL) is obtained ji^. To combine this result with data from oscillation experiments we use 
the x 2 from cosmology for four neutrinos, shown in Fig. 6 of Ref. |28| . 

In the following analysis we consider (3+1) schemes of the hierarchical type (first two 
spectra in Fig. P) and neglect the masses of the three lower mass states, such that m u ~ 



^ m LSND- 3 The A^cosmoiogy fr° m Ref- under this assumption is shown in the right 
panel of Fig. El as a function of AmJ SND . From the left panel of this figure one finds that 
cosmology excludes the region Am^ SND >3 eV 2 . However, the islands at Am 2 SND ~ 0.91 eV 2 
and 1.7 eV 2 are practically unaffected by cosmological data. Performing a GOF analysis 
shows that adding cosmological data gives Xpg = l^.l and 19.6 for the two islands. This 
should be evaluated for 5 DOf[I^, leading to a PG of 0.43% and 0.15%, respectively, 
which is comparable to the value given for (3+1) in Tab. ^ We note, however, that a 
modest improvement of the bound on m u by future cosmological data will drastically 
worsen the fit of the last remaining islands in the four-neutrino parameter space. 

A further cosmological constraint on four-neutrino schemes comes from Big Bang 
nucleo synthesis fBBN) |22l. l33j | . In the absence of any non-standard mechanism, like a 
large lepton number asymmetry, (3+1) as well as (2+2) schemes will lead to an effective 
number of neutrino species N u = 4]34]. This is in conflict with the upper bound on N v 

3 Note that all other four neutrino schemes will be much more disfavoured, since more mass states will 
contribute at the LSND mass scale. 
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Figure 5: Left: Allowed regions at 90% and 99% CL for (3+1) schemes without (solid 
and dashed lines) and including data from cosmology (coloured regions). The grey region 
is the 99% CL region of LSND0. Right: A% 2 from cosmological data|28| as a function 
of Aro? SND . 



derived in recent analyses [281. 35l . 36( including CMB and BBN data. Taking into account 
the uncertainty related to the primordial abundances of He and D, upper bounds in the 
range N v < 3.1 — 3.3 (2cr) 36j are obtained. 



7 Conclusions 

Performing a global analysis of current neutrino oscillation data we find that four-neutrino 
schemes do not provide a satisfactory fit to the data: 

• The strong rejection of non-active oscillation in the solar+KamLAND and atmo- 
spheric neutrino data rules out (2+2) schemes, independent of whether LSND is confirmed 
or not. Using an improved goodness of fit method especially sensitive to the combination 
of data sets we obtain a GOF of only 1.6 x 10~ 6 for (2+2) schemes. 

• (3+1) spectra are disfavoured by the disagreement of LSND with short-baseline 
disappearance data, leading to a marginal GOF of 5.6 x 10~ 3 . If LSND should be confirmed 
we need more data on v e and/or SBL disappearance to decide about the status of (3+1). 

• Recent cosmological data lead to further restrictions of the allowed parameter space. 
Only two islands aground AmJ SND ~ 0.91 eV 2 and 1.7 eV 2 in hierarchical (3+1) schemes 
survive. 

Our analysis brings the LSND hint to a more puzzling status, and the situation will 
become even more puzzling if LSND should be confirmed by the up-coming MiniBooNE 
experiment 0] . Recently an explanation in terms of five neutrinos has been proposed [38^. 
In such a (3+2) scheme solar and atmospheric data are explained dominantly by active 
oscillations, similar to the (3+1) case. The fit of SBL data is improved with respect to 
(3+1) by involving two large mass splittings Am^ ~ 0.9 eV 2 and Am^ ~ 20 eV 2 . Alter- 
native explanations of LSND by anomalous muon decayjl^] seem to be in disagreement 
with the negative result of KARMEN, or require a very radical modification of standard 
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physics, like CPT violation [ioj. 
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